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In any quantum or wave system dissipation leads 
to decoherence. Therefore, it was surprising in 
first instance when experiments on strongly lossy 
random lasers showed unambiguously by mea- 
surements of the photon statistics and of the las- 
ing mode volume that coherent feedback is pos- 
sible in such systems^'^'^. In coherent-feedback 
lasers the photons form a far-from-equilibrium 
condensate in the sense that a single quantum 
state is occupied by a macroscopic number of 
photons-*^^. We demonstrate that the lossy dy- 
namics of random lasers alone imply a finite lasing 
mode volume, thus resolving the puzzle about co- 
herent feedback without resonator. Our theory of 
random lasing including nonlinear gain and gain 
saturation predicts a characteristic dependence of 
this lasing mode volume on the pump intensity, 
which can be tested experimentally. The con- 
cept of a random laser (RL), i.e., lasing in a homoge- 
neously disordered, laseractive medium without cavity, 
was introduced early on by Letokhov. It was then fore- 
seen that light would be amplified by stimulated emission 
while diffusing through the random medium and cover- 
ing an unlimited region in space. Consequently, a con- 
tinuum of all the modes whose frequency lies within the 
laser transition line would participate in the lasing, lead- 
ing to an uncorrelated (Bose-Einstein) photon number 
statistics (distributed feedback) and homogeneous laser 
emission from the entire system volume. By contrast, it 
was found in experiments on ZnO powder lasers that the 
laser light emission occurs from spatially strongly con- 
fined lasing spots with sharp, discrete emission lines and 
Poissonian photon number statistics, all three features 
being unambiguous signatures of coherent feedback, i.e. 
of laser emission from a single, spatially confined mode^. 
This spatial confinement is at the heart of understanding 
coherent feedback and the formation of a non-equilibrium 
photon condensate state in homogeneous random lasers. 
However, its origin remained obscure ever since its discov- 
ery, because the commonly known, possible confinement 
mechanisms cannot explain the observed experimental 
facts. Anderson localization (AL) of light in disordered 
media can be ruled out, because the photonic scatter- 
ing mean free path £ is much longer than the wavelength 
of light A, so that AL does not occur in the systems 
at hand, and the loffe Regel criterion is not met, i.e. 
kl > 1. Random microcavities, preformed accidentially 
in the disordered medium, are unlikely to be the origin 



of coherent feedback, because neither would their intrin- 
sic surface roughness allow for localized cavity modes to 
exist nor could they explain the observed dependence of 
the average lasing spot radius on the pump intensitj^i^i^. 

Here we present a genuinely new mechanism for gener- 
ation of a finite coherence or mode volume and coherent 
feedback, based on dissipation at the system surface. We 
show that optical losses at the surface imply that the 
light field in the bulk of the system can be coherent (cor- 
related) only over a finite length ^. Hence, this defines a 
volume within coherent feedback occurs, leading to Pois- 
sonian photon statistics. We identify f with the average 
radius of a lasing mode. It is dynamically generated by 
the highly non-linear, lossy dynamics of open laser sys- 
tems and, as such, exists only in the lasing state. As 
we perceived, this confinement mechanism is closely re- 
lated to the causality of both the propagation of light 
and the transport of light intensity. Moreover, this being 
a dynamical effect, we predict ^ to have a characteristic, 
namely power law, dependence on the laser pump inten- 
sity, which is found to be in qualitative agreement with 
available experiments^ii^i^. 

As a RL system we consider specifically a layer of 
compressed powder of laser-active material, whereby the 
grains act as random scatterers and at the same time as 
amplifying medium. The layer has a thickness d and ex- 
tends infinitely in the x-y plane, as depicted in Fig. [T^. 
The pump light covers a wide surface area and affects the 
entire layer, so that the pump intensity may be assumed 
homogeneous across the entire volume of the considered 
system. The laser material is characterized by an atomic 
four-level scheme with transition rates 7^ as defined in 
Fig. Hb, although any other laser type would be possible. 

In order to address the problem of the coherence vol- 
ume in random lasers it is essential to set up the equa- 
tion of motion for the coherent part of the radiation. 
This part is gerenerated by stimulated emission only, al- 
though the quantum dynamics of the gain medium in- 
cludes also spontaneous emission. The time evolution 
of the electric radiation field can be parameterized by 
as E(t,r) — Eo(i, r) exp[—iujt], where uj/2tt is the light 
frequency and 'Eo{t,r)exp[—iuit] an amplitude function 
varying slowly on the scale of 27r/a;. We now observe that 
random lasers have to be treated as open systems, i.e., the 
loss-induced damping time is always much shorter than 
the lifetime 712 of the upper atomic laser level (class B 
laser). In this strongly damped regime, the stimulated 
part of the polarization follows the electric field instan- 
taneously. Hence, the wave equation for the stimulated 
emission is local in time and can be written in terms of 
the dielectric function £(r,Eo(i,r)) as 

[£i^:iMl^^^ + V^]E(t,r) = 0. 

Note that e(r, Eo(i,r)) incorporates the full, nonlinear 
laser dynamics through its dependence on the field am- 
plitude Eo(i, r) and depends on position r both explicitly 
because of the random position of the dielectric scatterers 
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and implicitly through its dependence on Eg. The las- 
ing state is characterized by a negative imaginary part 
of e(r, Eo(t, r)). It is a priori not known and therefore 
determined by solving a four level laser rate equation 
system (Siegmann) in steady state (see supplementary 
informations). 
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The Bethe-Salpeter-Equation is the approbiate equa- 
tion of motion for the coherent part of light intensity 
$. The Green's function G itself describes the behav- 
ior of the elctromagnetic wave, whereas the so called 
irreducible vertrex 7 reveals all interactions and the 
selfconsistence-loop is closed by the selfdependency of 
It has to be stressed that the propagation of light inten- 
sity in diffusive or strong scattering regimes obeys the 
second law of thermodynamics and thus the time scales 
of wave propagation and intensity propagation in fact 
can be separated. The results drawn from these consid- 
erations turned out to be really promising, because the 
diffusive processes, the interferences and localization pro- 
cesses can be clearly identified as separate contributions 
in the thereby derived diffusion constant D{^}) 



The electronic transition rates from the ground state 
to the uppermost state and from the lower lasing level 
to the ground state are assumed to be large compared 
to all other transitions. Thus we focus on the de- 
scription of the photon number Uph and the inversion 

Tin — No /Not — — T J'' I TTT according to the 

Einstein rate equations in steady state. 7p equals the 
pump rate, jnr represents the nonradiative decays and 
721 = 1/t2i features the stimulated emission rate, which 
is the inverse of the lifetime or relaxation time of the 
upper lasing level. The rsulting inversion of the occupa- 
tion number thus defines the microscopic optical gain and 
the challenge is now to relate this gain to the coherent 
(correlated) part of the diffusing radiation. Therefor the 
laser rate system is coupled to sophisticated model for 
microscopic transport of light in disordered amplifying 
random mediaF^"^ based on Vollhardt-Wolfle theory 
for the transport of electronsJ^, which allows for the con- 
sideration of the influences of the microstructure of the 
amlifying media. For the completeness of the descrip- 
tion of ramdom lasing it has been demanded earlier by 
Florescu et. al.^'^ that the effects due to multiple scat- 
tering of light and the internal Mie-resonances during 
these scattering events must be considered. We included 
those contributions by generalizing the diagrammatic de- 
scription for the transport of light and we particularly fo- 
cussed on the influence of interferences. Therefore we had 
to deal with a Dyson description of the propagation of 
the electromagnetic wave on the one hand and the Bethe- 
Salpeter-Equation for the propagation and dissipation of 
the energy density of light on the other. The complete 
description can be found in the supplementary, but to 
introduce the reader breefly to this fascinating subject, 
the Bethe-Salpeter-Equation for the two particle Green's 
function or in other words the energy propagator $ is 
given here - which looks a bit denaunting at first glance 
but actually it is a treasure chest for the random laser 
theory. 



D{n) [1 - i nReeujT^] = D^°* - T^D{n)M{uj). 

The first term on the right hand side _Dq°* = Dq + 
Db + Ds could be slightly overlooked on the search for 
interference effects which are represented through M{uj), 
the so called memory kernel. A detailed study reveals, 
however, that the memory kernel cannot cause true local- 
ization, i.e. I? = 0, in laser active media, and Dq°* cannot 
be compensated. Dq refiects the bare diffusion, the dis- 
sipative renormalizations due to absorbtion and/gain in 
the system are summarized in the term Df, and Dg. Af- 
ter all we remark that this diffusion constant for systems, 
which are in any way laser-active, will never vanish, and 
therefore AL of light in such systems seems to be beyond 
the question in the discussion of random lasing. 

From the microscopically derived correlation length 
^ (see supplementary information) in steady state for 
strong pumping (see Fig. 4) one can draw several sig- 
nificant conclusions for the characteristics of a random 
laser 



where a is a numerical constant an P the pumping. It 
is found that the natural antagonists, the pump strength 
P and the loss at the surfaces play a severe though in first 
instance counterintuitive role. The coherence length ^ in 
steady state clearly shrinks with the increase of the pump 
strength P whereas dissipation or loss at the surfaces 
of the depicted slab geometry (see Fig. 1) in contrary 
increases the correlation length. The meaning of both 
results become obvious when we go back to the roots of 
the idea of random lasing again. Dicke in 1968 proposed 
the possibility of the photon bomb by means of an infinite 
increase of the k-modes in an amplifying media. Our 
results in contrast proove that the nature of the onset of 
lasing in amplifying random media is clearly described by 



3 



a finite correlation length and thus a finite mode volume 
which marks the exponential decay of the spatially and 
spectrally coherent light intensity inside the ZnO powder 
slab (Fig 3b). Loss at the surface however leads to an 
increase of the correlation length ^. Both aspects, the 
dependency on the pumping and on the loss lead to gain 
saturation which govern the correlation volume in steady 
state. 

These results clearly forbid the development of a pho- 
ton bomb in lossy, finite, but amplifying random systems 
on the one hand, but also guarantee the coherence of the 



emission. The relation between the atomic population 
inversion and the dynamics of the coherent part of the 
light field described by the Bethc-Salpeter-Equation con- 
stitutes the parameter-free link between microscopic and 
mesoscopic processes in random lasing which has been 
missing up to the present. The description of growth 
of the photon number density in both approaches estab- 
lishes the selfconsistency chain that incorporates the full 
dynamics of the random laser. Hence the first parameter- 
free description of random lasing is established. 
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FIG. 2; The full and self-consistent diffusion constant D in 
units of the bare diffusion constant Do is shown as a function 
of the width Z across the sample for pump rate P/721 = 2. 
The sample width is set to be d = 40ro and d = 80ro respec- 
tively as indicated in the legend , where ro is the radius of the 
scatterers. As discussed in the text, the ratio D/ Do remains 
practically unchanged, since the relative change observed in 
the graph is of the same order as the accuracy, involved in the 
numerical evaluation. This implies that interference effects do 
not change as a function of the sample width. The diffusion 
constant D can be seen to behave inversed compared to the 
inversion 712. Both quantities show a weak dependency to the 
sample's width. 
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FIG. 3: a The correlation volume of the RL mode strongly 
decays with encreasing pump strength, b Calculated photon 
density 




FIG. 4: Correlation length ^ as a function of the inverse 
square root of the pump rate, measured in units of the in- 
verse square root of the transition rate 721, at the sample 
surface, i.e. Z = 0.5d. 

The blue lines serve as a guide to the eye, emphasizing the lin- 
ear behavior of the correlation length in this plot. This clearly 
reveals an inverse square root behavior of ^ as a function of 
pumping above threshold. 



